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The invariant mass spectrum of the lepton pair in inclusive semilep- 
tonic B — » XJlV decay yields a model independent determination of 
\V u b\ [[![. Unlike the lepton energy and hadronic invariant mass spec- 
tra, nonperturbative effects are only important in the resonance region, 
and play a parametrically suppressed role when dT/dq 2 is integrated over 
q 2 > (m_B — mo) 2 '', which is required to eliminate the B — > X c £v back- 
ground. We discuss these backgrounds for q 2 slightly below {m B — m D ) 2 , 
and point out that instead of q 2 > (m B — m D ) 2 = 11.6 GeV 2 , the cut 
can be lowered to q 2 > 10.5 GeV 2 . This is important experimentally, 
particularly when effects of a finite neutrino reconstruction resolution are 
included. 
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A precise and model independent determination of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element V u b is important for testing the Standard Model at B factories 
via the comparison of the angles and the sides of the unitarity triangle. 

If it were not for the huge background from decays to charm, it would be straight- 
forward to determine \V u b\ from inclusive semileptonic decays. Inclusive B decay 
rates can be computed model independently in a series in and a s (nib) using 

an operator product expansion (OPE) [|2],|3],|4|,|5|] , and the result may schematically be 
written as 

dr= p q uarkX r +A + /(A 1 ^) + ___ + a £ a| y 
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At leading order, the B meson decay rate is equal to the b quark decay rate. The 
leading nonperturbative corrections of order Ag^/m 2 . are characterized by two heavy 
quark effective theory (HQET) matrix elements, usually called Ai and A2. These 
matrix elements also occur in the expansion of the B and B* masses in powers of 
Aqcd/W 

m B = m b + A \- ... , m B * = m h + A V ... ■ (2) 

Znib 2nib 

Similar formulae hold for the D and D* masses. The parameters A and Ai are inde- 
pendent of the heavy b quark mass, while there is a weak logarithmic scale dependence 
in A2. The measured B* — B mass splitting fixes ^(mb) = 0.12 GeV 2 , while A and 
Ai (or, equivalently, a short distance b quark mass and Ai) may be determined from 
other physical quantities Thus, a measurement of the total B — > X u iT7 rate 

would provide a ~ 5% determination of \V u b\ P,|10"l|. 

Unfortunately, the B — > X u tP rate can only be measured imposing cuts on the 
phase space to eliminate the ~ 100 times larger B — > X c tP background. Since the 
predictions of the OPE are only model independent for sufficiently inclusive observ- 
ables, these cuts can destroy the convergence of the expansion. This is the case for 
two kinematic regions for which the charm background is absent and which have re- 
ceived much attention: the large lepton energy region, Et > (m B — m 2 D )/2niB, and 



the small hadronic invariant mass region, nix < m D [JTT],[T2],[T3|] . 

The poor behaviour of the OPE for these quantities is slightly subtle, because 
in both cases there is sufficient phase space for many different resonances to be pro- 
duced in the final state, so an inclusive description of the decays is still appropriate. 
However, in both of these regions of phase space the B — > X u tP decay products are 
dominated by high energy, low invariant mass hadronic states, 

E x ~ m b , m 2 x ~ A QC D"if> > Aq CD (3) 

(where Ex and nix are the energy and invariant mass of the final hadronic state). In 
this region the differential rate is very sensitive to the details of the wave function of 
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the b quark in the B meson. Since the OPE is just sensitive to local matrix elements 
corresponding to expectation values of operators in the meson, the first few orders in 
the OPE do not contain enough information to describe the decay, and as a result 
the OPE does not converge. 

This is simple to see by considering the kinematics. A b quark in a B meson has 
momentum 

p» = m bV » + p (4) 

where is the four-velocity of the quark, and k^ is a small residual momentum of 
order Aqcd- If the hadron decays to leptons with momentum q and light hadrons 
with total momentum px, the invariant mass of the light hadrons may be written 

m\ = (m b v + k -q) 2 = (m b v - q) 2 + 2k ■ (m b v - q) + 0(Aq CD ). (5) 

The first term in the expansion is 0(m 2 ) over most of phase space, while the second 
is 0(Aqcd"^), and so is suppressed over most of phase space. The OPE presumes 
that this power counting holds, so that the second term may be treated as a small 
perturbation. However, if Ex is large and mx is small, m b v — q is almost light-like, 

rrw? - q» = (E x , 0, 0, E x ) + 0(A QCD ) (6) 

in the b rest frame where t> M = (1,0,0,0). Since Ex ~ 0(m b ), (m b v — q) 2 = 
O(Aqcd^)- Thus, in this region the first two terms in (f|) are of the same order 
(but still parametrically larger than the remaining terms), and the invariant mass of 
the final hadronic state reflects the distribution of the light-cone component of the 
residual momentum of the heavy quark in the hadron, 

m 2 x = (m b v - q) 2 + 2E x k + + k + = k + k 3 . (7) 

Since the differential rate in this region depends on the invariant mass of the final 
state, it is therefore sensitive at leading order to the light-cone wave function of the 
heavy quark in the meson, f(k+). 

In terms of the OPE, this light-cone wave function arises because of subleading 
terms in the OPE proportional to ExA.qcT)/m x , which are suppressed over most 
of phase space but are 0(1) in the region @. It has been shown that the most 
singular terms in the OPE may be resummed into a nonlocal operator whose matrix 
element in a B meson is the light-cone structure function of the meson. Since f(k + ) 
is a nonperturbative function, it cannot be calculated analytically, so the rate in the 
region (|3|) is mo del- dependent even at leading order in Aq CD /m b . 

The situation is illustrated in Fig. ffl, where we have plotted the lepton energy and 
hadronic invariant mass spectra in the parton model (dashed curves) and incorporat- 



ing a simple one-parameter model for the distribution function (solid curves) JT7] 
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Figure 1: The shapes of the lepton energy and hadronic invariant mass spectra. The 
dashed curves are the b quark decay results to 0(a s ), while the solid curves are obtained 
by smearing with the model distribution function f(k + ) in Eq. (||). The unshaded side of 
the vertical lines indicate the region free from charm background. 



The differences between the curves in the regions of interest indicate the sensitivity of 
the spectrum to the precise form of }'{k + ). Currently, there are measurements of \V U },\ 
from both methods. From the lepton energy cut, the PDG reports (V^/V^l = 0.08 ± 
0.02, while a recent DELPHI measurement using the hadronic invariant mass cut 
gives \ V ub /V cb \ = 0.1031^ (syst.) ± 0.016 (stat.) ±0.010 (theory) p|. In both cases, 
the theoretical error is an estimate based on varying different models of f{k + ), and so 
these measurements are no more model-independent than the exclusive measurement 
from B — > ptv. While it may be possible in the future to extract f{k + ) from the 
B — > X s 7 photon spectrum |14|,|T8|, unknown order Aqcd / m b corrections arise when 



relating this to semileptonic b — > u decay, limiting the accuracy with which \ V u b\ may 
be obtained. 

Clearly, one would like to be able to find a cut which eliminates the charm back- 
ground but does not destroy the convergence of the OPE, so that the distribution 
function f{k + ) is not required. In Ref. |L] we pointed out that this is the situation for 
a cut on the dilepton invariant mass. Decays with q 2 = (pi+Pv) 2 > {m B —m D ) 2 must 
arise from b u transition. Such a cut forbids the hadronic final state from moving 
fast in the B rest frame, and simultaneously imposes mj < and Ex < m^. Thus, 
the light-cone expansion which gives rise to the shape function is not relevant in this 
region of phase space |T3]JT9|1 . The effect of smearing the q 2 spectrum with the model 
distribution function in Eq. (g) is illustrated in Fig. Q It is clearly a subleading effect. 
The Dalitz plots relevant for the charged lepton energy and hadronic invariant mass 
cuts are shown in Fig. Note that the region selected by a q 2 cut is entirely contained 
within the m 2 x cut, but because the dangerous region of high energy, low invariant 
mass final states is not included with the q 2 cut, the OPE does not break down. It is 
also important to note, however, that the q 2 cut does make the OPE worse than for 
the full rate; as we will show, the relative size of the unknown Aq CD /t?t| terms grows 
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Figure 2: The dilepton invariant mass spectrum. The notation is the same as in Fig. |l]. 



the effective expansion 



as the q 2 cut is raised. Equivalently, as was stressed in 
parameter for this region is Aq CD /m c , not AQ CD /m b . 

The B — > X u lv decay rate with lepton invariant mass above a given cutoff can 
therefore be reliably computed working to a fixed order in the OPE (i.e., ignoring the 
light-cone distribution function), 
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where q 2 = q 2 /m 2 , /?o = 11 — 2nj/3, and T 
b — > u decay rate. The ellipses in Eq. 

a: 
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F \ V u b\ m|/(192 7r) is the tree level 
denote terms of order (Aqcd/^) 3 and order 
terms not enhanced by (3q. The function X(q 2 ) is known analytically |HJ], whereas 



Y(q 2 ) was computed numerically [^2]. The order 1/mf nonperturbative corrections 
are also known as are the leading logarithmic perturbative corrections propor- 
tional to a™ log n (m c /mb) [pQ[| . The matrix element of the kinetic energy operator, Ai, 
only enters the q 2 spectrum in a very simple form, because the unit operator and the 



kinetic energy operator are related by reparameterization invariance [23 . 

The relation between the total B — > X u lv decay rate and \V ub \ is known at the 
~ 5% level @,[TQ|, 



\V ub \ = (3.04 ± 0.06 ± 0.08) x 10" 



'B{B -> x u 



1.6 ps 
0.001 x F{q 2 ) 



1/2 



(10) 



where F(q 2 ) is the fraction of B — > X u iv events with q 2 > q 2 , satisfying -F(O) 



1. The errors explicitly shown in Eq. (|10|) are the estimates of the perturbative 
and nonperturbative uncertainties in the upsilon expansion |§ respectively. At the 
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Figure 3: Dalitz plots relevant for B —* X u lV^. The shaded regions indicate the part 
of phase space where B — » X c tp£ background is present, and the vertical dashed regions 
corresponds to the cut q 2 > (mg — mo) 2 - In the q 2 — qo plane, the horizontal dashed region 
corresponds to an invariant mass cut m\ > m 2 D , whereas in the q 2 — Eg, plane the horizontal 
dashed region corresponds to the charged lepton energy cut Eg > {m 2 B — m D )/2m,B- Note 
that at tree level, b — > u semileptonic decay populates the entire triangle on the right-hand 
plot, but only the right boundary of the left-hand plot. 



present time the biggest uncertainty is due to the error of a short distance b quark 
mass, whichever way it is denned [p0[| . (This can be cast into an uncertainty in an 
appropriately denned A, or the nonperturbative contribution to the T(IS') mass, etc.) 
By the time the q 2 spectrum in B — > X u £v is measured, this uncertainty should be 
reduced from extracting nib from the hadron mass |J or lepton energy spectra in 
B — ► X C £V, or from the photon energy spectrum |§ in B — > X s ^. The uncertainty 
in the perturbation theory calculation will be largely reduced by computing the full 
order a 2 correction in Eq. (|i~0|) . The largest "irreducible" uncertainty is from order 
AQ CD /mjj terms in the OPE, the estimated size of which is shown in Fig. f|, together 
with our central value for F(q^), as functions of qfi. 

There is another advantage of the q 2 spectrum over the nix spectrum to measure 
\V u b\. In the variable mx, about 20% of the charm background is located right next 
to the b — > u "signal region", mx < mD, namely B — > Div at nix — m D- I n the 
variable q 2 , the charm background just below q 2 = (mg — ni£>) 2 comes from the 
lowest mass X c states. Their q 2 distributions are well understood based on heavy 
quark symmetry | 25| , since this region corresponds to near zero recoil. Fig. |5] shows 
the B — > DiV and B — > D*iV decay rates using the measured form factors (and 
\V u b\ = 0.0035). The B — > X u £u rate is the flat curve. Integrated over the region 
q 2 > (niB — ma*) 2 ~ 10.7 GeV 2 , the uncertainty of the B — > D background is small 
due to its (w 2 — 1) 3//2 suppression compared to the B — > X u £v signal. This uncertainty 
will be further reduced in the near future. This increases the b — > u region relevant 
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Figure 4: (a) The fraction of B — > X U £U events with q 2 > F^q^), in the upsilon expansion. 
The dashed line indicates the lower cut q$ = {tub — mi)) 2 ~ 11.6 GeV 2 , which corresponds 
to F = 0.178 ± 0.012. The shaded region is the estimated uncertainty due to Aq CD /m^ 
terms; which is shown in (b) as a percentage of F^q^). 



<K 

dq 2 



B^Dlv 



11 11.2 11.4 

q 2 (GeV 2 ) 



B^XJV 



Figure 5: Charm backgrounds near q 2 = (tub — mo) 2 (arbitrary units). The shaded region 
denotes the uncertainty on the B — > Dtu rate. 



for measuring \ V u i\ by ~ 1 GeV 2 . The B — ► D* rate is only suppressed by (w 2 — l) 1 / 2 
near zero recoil, and therefore it is more difficult to subtract it reliably from the b — > u 
signal. The nonresonant Dir final state contributes in the same region as B — > D*, 
and it is reliably predicted to be small near maximal q 2 (zero recoil) based on chiral 
perturbation theory |27|]. The D** states only contribute for q 2 < 9 GeV 2 , and some 
aspects of their q 2 spectra are also known model independently . 

Concerning experimental considerations, measuring the q 2 spectrum requires re- 
construction of the neutrino four- momentum, just like measuring the hadronic invari- 
ant mass spectrum. A lepton energy cut may be required for this technique, however, 
the constraint q 2 > (tjib — mc) 2 automatically implies Eg > (thb — mo) 2 /^tub — 
1.1 GeV in the B rest frame. Even if the Eg_ cut has to be slightly larger than this, 
the utility of our method will not be affected, but a calculation including the effects 
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of arbitrary En and q 2 cuts would be required. If experimental resolution on the 
reconstruction of the neutrino momentum necessitates a significantly larger cut than 
q 2 = (itlb — ttld) 2 , then the uncertainties in the OPE calculation of F(q^) increase. 
In this case, it may be possible to obtain useful model independent information on 
the q 2 spectrum in the region q 2 > mLm ~ 13.6 GeV 2 from the q 2 spectrum in the 
rare decay B — > X s £ + £~ , which may be measured in the upcoming Tevatron Run-II. 

In conclusion, we have shown that the q 2 spectrum in inclusive semileptonic B — > 
X u lu decay gives a model independent determination of \V u b\ with small theoretical 
uncertainty. Nonperturbative effects are only important in the resonance region, and 
play a parametrically suppressed role when dT/dq 2 is integrated over q 2 > (m# — 
mo) 2 , which is required to eliminate the charm background. This is a qualitatively 
better situation than other extractions of \V u b\ from inclusive charmless semileptonic 
B decay. 
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